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1. Introduction 

It has been known for many years that the Hamilton-Jacobi equation of classical me- 
chanics defines a function which appears to be the eikonal of a wave equation, and therefore 
that classical mechanics appears to be a ray approximation to some wave theory 1 . The 
propagation of rays of waves in inhomogeneous media appears, from this point of view 
(as a result of the application of Fermat's principle), to correspond to geodesic motion 
in a metric derived from the properties of the medium 2 . This geometrical interpretation 
has been exploited recently by several authors to construct models which exhibit three 
dimensional analogs of general relativity by studying the wave equations of light in an 
inhomogeneous medium 3 , and, to achieve four dimensional analogs, sound waves and elec- 
tromagnetic propagation in inhomogeneously moving materials 4 . Visser et a/ 5 ' 6 ' 7 have 
pointed out that condensed matter systems such as acoustics in flowing fluids, light in 
moving dielectrics, and quasiparticles in a moving superfluid can be used to mimic kine- 
matical aspects of general relativity. Leonhardt and Piwnicki 8 and Lorenci and Klippert 9 , 
for example, have discussed the case of electromagnetic propagation in moving media. In 
order to achieve four dimensional geodesic flows, it has been necessary to introduce a mo- 
tion of the medium. * There is considerable interest in extending these analog models for 
the kinematical aspects of gravity to include dynamical aspects, i.e., considering gravity 
as an emergent phenomenon 5 ' 6 . 

The manifestly covariant classical and quantum mechanics introduced by 
Stueckelberg 12 in 1941 has the structure of Hamiltonian dynamics with the Euclidean 
three dimensional space replaced by four-dimensional Minkowski space (since all four of the 
components of energy-momentum are kinematically independent, the theory is intrinsically 
"off-shell"). The dynamical evolution of the system is governed by a "world time" r. This 
theory leads to five dimensional wave equations for the associated gauge fields. The self- 
interaction problem of the relativistic charged particle has been recently studied 13 , where 
it was shown that the radiation field is associated with excursions from the mass shell; 
highly non-linear terms appear in the resulting generalized Lorentz force. In the limit that 
the motion maintains a state very close to the mass shell limit, the equation reduces to 
that of Dirac 14 . 

The known bound state spectra and corresponding wave functions for the (spinless) 
two-body problem in potential theories formulated in a manifestly covariant way have been 
worked out 15 . The classical relativistic Kepler orbits have been studied in detail 12 ' 16 . 

In order for the Stueckelberg-Schrodinger equation of the quantum form of this theory 
to be gauge invariant, it is necessary to introduce a fifth gauge field, compensating for the 
derivative with respect to the invariant evolution parameter r 17 > 18 . Generalized gauge 
invariant field strengths, with q,p = 0, 1, 2, 3, 5 occurring in the Lagrangian to second 
order generate second order field equations analogous to the usual Maxwell equations, 



* We remark that Obukhov and Hehl 10 have shown that a conformal class of metrics 
for spacetime can be derived by imposing constrained linear constitutive relations between 
the electromagnetic fields (E, B) and the excitations (D, H), using Urbantke's formulas 11 , 
developed to define locally integrable parallel transport orbits in Yang-Mills theories (on 
tangent 2-plane elements on which the Yang-Mills curvature vanishes). 
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with source given by the four-vector matter field current and an additional Lorentz scalar 
density. The canonical second quantization of this theory was studied in ref. 19. Taking 
the Fourier transform of these equations over the invariant parameter, as we demonstrate 
below, one seees that the zero frequency component (zero mode) of the equations coincides 
with the standard Maxwell theory (the fifth field decouples). The Maxwell theory is 
therefore properly contained in the five-dimensional generalization as the zero frequency 
component. The general form of the theory has variations as a function of world time, 
about the (world time independent) Maxwell theory. 

In the quantum case the four-currents are given by bilinears in the wave function 
containing first derivatives, and the fifth source is the (Lorentz) scalar probability density. 
The symmetry of the homogeneous equations, which can be 0(3, 2) or 0(4, 1), depending 
on the sign chosen for raising and lowering the fifth index; it is not realized in the inhomo- 
geneous equations, since the spacetime current components and the fifth component (the 
quantum mechanical probability density \i^>\ 2 ) cannot be transformed into each other by 
a linear transformation. Hence, without augmenting the symmetry of the matter fields 
beyond 0(3, 1), the fifth field, whose source in the Maxwell-like equations is the probabil- 
ity density (or, classically,the matter density) in spacetime, can play a special role. There 
appears to be no kinematic basis for choosing one or the other of these signatures; atomic 
radiative decay, for example, contains points in phase space (for radiation of off-shell pho- 
tons) for either type. We note, however, that the homogeneous equations corresponding to 
the 0(4, 1) signature appear, under Fouries transform of the r variable, as Klein-Gordon 
type wave equations with positive mass-squared (physical particles) , but for the 0(3,2) 
choice of signature, these equations have the "wrong" sign (tachyonic) for interpreting the 
additive term as a mass-squared. As a physical example of how both metrics may play a 
role, let us suppose that the off-shell radiation impinging on an atomic nucleus, for exam- 
ple, has the 0(3, 2) signature. During the interaction, as the state evolves as a function of 
r, the metric may make a transition to the 0(4, 1) form. This transition could correspond 
to the observed phenomenon of photoproduction (the vector dominance model). We there- 
fore leave open the question of a definitive choice of the signature for the five dimensional 
radiation field at this stage. 

In either case, the four dimensional spacetime submanifold coincides with that of 
Minkowski; on this manifold, the fields are defined as a family of functions paramaterized 
by r. This parameter controls the evolution generated by the Hamiltonian of the system 
of particles and radiation. That this parameter occurs in the wave equations with timelike 
or spacelike metric does not change its physical interpretation. At each value of r, the four 
dimensional configuration corresponds to a Maxwell wave; these waves, however, evolve 
non-trivially in r, reflecting the dynamics of the system. For example, in classical charged 
particle scattering, the conserved currents are constructed (as we remark below) in the 
four-dimensional Maxwell theory, by integration of a local current j^S (x — x{s)) over all 
proper time s (which we may identify here with the world time). The trajectory x fl (s) is 
not, however, known until the problem is solved; the system is not a priori integrable. In 
the higher dimensional theory we discuss here, equations of motion determine the trajectory 
x(t), and in the neighborhood of a given r, the generalized (sometimes called pre-Maxwell) 
equations determine the fields as functions of r. This system is explicitly, in principle, 
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integrable. The integral of the results over all r (on both fields and currents) then provides 
exact solutions of the corresponding Maxwell problem, with currents computed from the 
resulting trajectories in the standard way. Since, however, the Lorentz force is non-linear, 
it is clear that the dynamics forming the trajectories may be very different from that of 
the usual Lorentz force 13 . 

The structure of the gauge theory obtained from the non-relativistic Schrodinger equa- 
tion is precisely analogous. The four dimensional gauge invariant field strengths obey 
second order equations, for which the sources are the vector currents and the scalar prob- 
ability density (or, classically, the matter density in three dimensional space). No linear 
transformation can connect the Schrodinger probability density with the vector currents, 
so that the scalar density and the fourth gauge field can play a special role. 

Since the eikonal approximation naturally lowers the dimension of the differential 
equations describing the fields by one, the eikonal approximation to the five-dimensional 
field equations results in four dimensional differential equations. In the presence of a non- 
trivial dielectric structure of the medium, the four dimensional field equations resulting 
from the eikonal approximation can describe geodesic motion in four dimensional space- 
time without the necessity of adding motion to the medium. We emphasize that the 
underlying manifold, on which the fields are defined, is a flat Cartesian space, but that 
the dynamically induced trajectories are curved, and can be described by the geodesies of 
a pseudo-Riemannian manifold. This result forms our basic motivation for studying the 
generalized dynamics of Stueckelberg 12 in this context. 

We start with a review of the simpler case of Maxwell wave propagation in a three 
dimensional inhomogeneous, anisotropic medium in which we introduce a method which 
permits us to identify clearly the structure of the effective mechanics which emerges in 
the eikonal approximation. In Section 3, we study the eikonal structure of waves in a 
5D inhomogeneous medium, in which Minkowski spacetime is embedded, and in Section 
4, a wave equation of Schrodinger type, and show that the resulting rays have a precise 
analog with the results of Kline and Kay 2 (who have studied the three dimensional ray 
limit of the usual Maxwell equations), but in a 4D spacetime manifold with a pseudo- 
Riemannian structure. Kline and Kay 2 show that the rays are geodesic in the metric 
associated with the anisotropic inhomogeneous medium. As for the case of Kline and 
Kay 2 , it follows from the existence of a Hamiltonian that the corresponding Lagrangian 
obeys an extremum condition which describes the rays as geodesies. We show that there 
is mass flow along these rays, and that the flow is controlled by generating functions of 
Hamiltonian type, establishing a relation between geodesic flow and a relativistic particle 
mechanics of symplectic form. 

2. Eikonals of the Maxwell Wave Equation 

We first consider the case of Maxwell wave propagation in a three dimensional inho- 
mogeneous, anisotropic dielectric media, as in ref. 2. We find the Fresnel surfaces in terms 
of quadratic forms which are the solutions of an eigenvalue equation, where each polariza- 
tion is associated with a Riemannian metric 2,20 (of Finsler type) giving rise to geodesic 
flow. The study of the eigenvalue equations, rather than the determinant condition used 
in ref. 2, permits us to identify more clearly the structure of the effective mechanics which 
emerges. 
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We start by writing Maxwell's equations in an inhomogeneous, anisotropic medium. 
These are of the usual form 

V-D = p (la) 
VxH- — =J (16) 
V • B = (lc) 

VxE + — = 0, (Id) 

where p is the charge density, J the current density of the sources, and, with i,j = 1, 2, 3 
(we assume here that the coordinates labelled by these indices are Cartesian, and we write 
indices up or down for convenience of expression), where, however, the relations 

Di = e tJ (x.)Ej , Bi = fiij(pc)Hj (2) 

reflect the properties of the medium (we assume all medium properties to be independent 
of time and that, for our present purposes, there is no mixing between electric and magnetic 
fields through the constitutive tensor relations). We then multiply Eq.(ld) with the matrix 
act with the curl operator and substitute Eq.(lb). We obtain, in the absence of 
sources, 

Vx(^(VxE)) + 5 = (3) 
Substituting for D using Eq.(2), one obtains 

e- 1 (Vx(^- 1 (VxE))) + ^=0 (4) 

We rewrite Eq. (4) in index notation (we use Latin indices for space components 1,2,3, 
and t for the time index): 

(e- 1 ) rl s lmn d m (^- 1 ) ni (e ijk d j E k ) + d 2 t E r = (5) 

Assuming a solution of the form E = Aexpiuj(t — y?(x)) and using the eikonal approxi- 
mation (for large u>), one obtains 

-{e- 1 ) rl e lmn {ii- 1 ) ni e l3k d m yd J yA k = A r . (6) 

For a given wave front at a given time (p(x.)\t=t = c o one can think of curves, origi- 
nating from each point on the wave front, and everywhere in the direction of the gradient 
of <p(x). These lines can be considered to be trajectories on which the different elements 
of the wave-front propagate. The trajectories of the wave-front normals, in an anisotropic 
medium, however, are in general not in the same direction at a given point as the direction 
of propagation of the energy of radiation, i.e., of the Poynting vector. We shall refer to the 
trajectories which follow the Poynting vectors as rays. We now show that the equation for 
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the rays is a geodesic equation. The metric determining the geodesies is fixed by the elec- 
tric and magnetic properties of the medium and, in general, depends on the polarization 
of the field. 

We denote the wave front gradient (which can be interpreted, as will be seen later, as 
the momentum flowing along the ray) 

Pm = d m lf. 

Equation (6) then takes the following form: 

-(e- 1 ) rl e lmn (ix- 1 ) ni e ijk p m p i A k = A r , (7) 

Multiplying Eq. (7) by p s e sr , the left hand side vanishes, and we obtain 

P ■ (eA) = 0, (8) 

the eikonal form of Eq.(la). The eikonal solution is therefore consistent with all of Maxwell 
equations in the non-homogeneous medium. We now define (obviously positive definite for 
e and fi scalar; the symmetric in mj part is symmetric in rk if e and fi are symmetric) 

M rk j = -(e- 1 ) rl e lmn (iJ,- 1 ) ni e ijk (9) 

Due to the relation (8), the 3x3 matrix M[ k piPj acts in a two dimensional subspace; 
one may replace A by v = eA if is replaced by e^M^e"^., which takes the two 

dimensional subspace orthogonal to p into itself. Eq.(7) may then be written as 

Ml k p l p 3 v k = v r , (10) 

where M = eMe -1 . Eq.(10) clearly imposes a condition on the magnitude of the vector 
p. As we shall see, it is restricted to two discrete values. 

If we express the tensor M^j i n coordinates for which the 3 direction is parallel to 
the momentum p, we obtain the eigenvalue condition 

Mgtfv" = v>r, (11) 

where the primed quantities are in the momentum oriented coordinate system; since v 1 
provides support only in the two dimensional subspace orthogonal to p, the matrix M^ q 
is nonzero only in that subspace. It is not, in general, degenerate. If the eigenvalues 
are A( a ), a = 1,2, it follows from (11) that p'^ = |p| 2 must have the values Ip^l 2 = 
A^. Transforming back to the original frame, we see that (10) can be satisfied only 
for two polarization modes v k a y and corresponding momenta of the same direction with 

magnitudes |p( a )| 2 = A^. Multiplying, for each a, both sides by (normalized) and 
summing on r, we obtain the form 

HM=pMpMgl? -1 = 0, (12) 
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where 

9^=V(a)K k v k {a) . (13) 

We shall identify p^ below as a canonical momentum, orthogonal to the surface defined 
by <p; the (p^ a ^ direction dependent) matrices gf^ therefore act as (Finsler type) metrics 
for each of the polarizations 20 . 

Since depends only on x and the direction of p^ at any given point in space, 
the condition (12) determines the magnitude of p^ a \ and therefore describes a surface. 
These surfaces, described for a = 1,2, coincide with the Fresnel surfaces defined by Kline 
and Kay 2 by means of the determinant of coefficients of the eigenvalue equation. We have 
examined here directly the eigenvalue equations since we shall follow this method in the 
5D case. It is shown for the three dimensional case in ref. 2 (the proof is given below 
for the similar four dimensional problem) that d ( a )H^ is parallel to the Poynting vector 

Pi 

(clearly the same direction for each a). If we parametrize the flow along a given ray at x 
with some parameter s, this statement can be written as 

i = <W_ = xd {a) 

ds p\ v ' 

where A is a scale on s. The total derivative of with respect to Xi is given by 

dH^ _ dH^ dH^ dp ( k a) 
dx i ~ dx l + Q p (a) dx i ^ ' 

This quantity must vanish, since the derivative relates neighboring Fresnel surfaces, on 
which (in this mode) is zero. Substituting (14) in (15) and using 

d Pk a) _ d 2( P _ d Pr a) (16) 



we obtain 



which gives 



dx l dx % dx k dx k ' 



dx l dx k 



p 4 (o) = -\d Xi H^. (18) 



Eqs. (14) and (18) correspond to the locally symplectic structure of a Hamiltonian flow 
generated by the function (12) in each mode. Moreover, one sees that the geodesic equation 
associated with the metric equivalent to this Hamiltonian flow. This result agrees 

with application of the Fermat principle. 



3. Eikonals of the 5D Wave Equation 

In this section we apply a technique similar to that used above to study the structure of 
wave equations in five dimensions which follow as a consequence of the requirement of gauge 
invariance of the covariant classical and quantum mechanics of Stueckelberg 12 ' 17-19 ' 21 . We 
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show that there is a Hamiltonian form for the generation of the rays in this case as well, and 
that these rays form spacetime geodesies in a metric space determined by the "dielectric" 
constitutive tensor relating the tensor fields (analogous to (E, B)) and their corresponding 
excitation fields (analogous to (D,H)). 

We shall use Greek letters for space time indices (/j = 0,1,2,3) and Latin letters 
to include a fifth index representing the Poincare invariant r parameter in addition to 
the usual 4 spacetime coordinates (e.g., q = 0,1,2,3,5). The analysis proceeds by a 
generalization of the method discussed in Section 2 above. The generalized electromagnetic 
field tensor is written 

fqi <?2 = &q 1 Clq 2 dq 2 Cl qi , 

where a q are the so-called pre-Maxwell electromagnetic potentials ( the fifth gauge potential 
as is required for gauge compensation of id$, generating the evolution of the Stueckelberg 
wave function 12 ' 17 ). 

We introduce the dual (third rank) tensor 

uh hh — c l i h h qi <?2 f 

c J qi <?2 1 

where e ll 12 13 91 q2 is the antisymmetric fifth rank Levi-Civita tensor density. The homoge- 
neous pre-Maxwell equations are then given by 

d h k hhh = 0, (19) 

or, more explicitly (the 5 index is raised with signature ±, according to, as discussed above, 
(9(4, 1) or 0(3,2) symmetry of the homogeneous field equations), 

d 5 e hh5q ^ 2 f qiq2 +d a e llh(7 ^ q2 f qiq2 =0. (20) 

We now divide Eq.(20) into two cases. In the first, we take the indices l±, I2 to correspond 
only to space-time indices: 

d 5 e^ 5X<7 h a + 2cU^ A5 /a5 = 0, 

or 

d 5 e^^h a + 2d a e^ x f X5 = 0, (21) 

where e fIuaX is the four dimensional Levi-Civita tensor density. This equation, on the 
0-mode (r independent Fourier components) does not involve any of the usual Maxwell 
fields but only the fifth (Lorentz scalar) electromagnetic field. The second set from Eq.(20) 
corresponds to /ior/2 = 5. It is clear then that all the other 4-remaining indices must be 
space-time indices and we obtain 

e^ 5v d a f 5v = - e^dafsu = 0. (22) 

It is this equation that reduces on integration, over all r, to the two usual homogeneous 
Maxwell equations. On the zero mode the r derivative disappears and therefore Eq.(21) 
reads for ji = , v = i, V x f = and for the \i = i , v = j components, d t f + V/o = where 
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we have called = (/o,f). Equation (22) then reduces to the standard homogeneous 
Maxwell equations; therefore we see that the homogeneous pre-Maxwell equations (i.e., 
Eq.(19)) for the 5 component fields and the spacetime component fields decouple on the 0- 
mode. This has the effect of reducing the pre-Maxwell system of equations, as we discuss 
below as well for the inhomogeneous equations, to the usual Maxwell equations. With 
appropriate identification of the integrated quantities, the zero mode of the pre-Maxwell 
equations coincides with the Maxwell theory (it is for this reason that the five dimensional 
gauge fields associated with the Stueckelberg theory are called "pre-Maxwell" fields). 

We now turn to the current-dependent pre-Maxwell equations. These can be written: 

d h n l2h =-j h , (23) 

where n^ 2 are the matter induced (excitation) fields (corresponding to H, D in the 4D 
theory). We remark that, restricting our attention to the spacetime components of Eq. 
(23), which then reads 

d b n& + d v n>"' = -f, (23') 

we may extract the 0-mode by integrating over all r. Since j k satisfies the five dimensional 
conservation law dkj k = 0, its integral over r (assuming 12 ' 17 f ~ * for r — > ±oo) reduces 
to the four dimensional conservation law <9 M J M = 0, where J M = J j M (x, r)dr is the 0-mode 
part of (this formula for the conserved is given in Jackson 22 ). The first term of the 
left side of (23') vanishes under integration (assuming that n M5 — > for r — > ±oo), and 
one obtains the form 

d v F^ = J", 

where we may identify the zero mode fields F^ u , as above, with the Maxwell fields H, D. 
With the zero mode of (22), we see that the Maxwell theory is properly contained in the 
five dimensional generalization we are studying here. We see that in Eq. (23) the I2 = \x 
component reduces on the 0-mode to the standard Maxwell equations and for I2 = 5 we 
get d t n° + V • n = (using similar notation as above, n M5 = (no,n)); thus for the pre- 
Maxwell equation involving the currents, on the 0-mode, the fields associated with the fifth 
component and the spacetime fields decouple. 

We now assume the existence of linear constitutive equations in the dynamical struc- 
ture of the 5D fields in a medium which connects the n tensor-field to the k tensor-fields us- 
ing a fifth rank tensor £ which is a generalization of the fourth rank covariant permeability- 
dielectric tensor 23 which relates the E, B fields to the excitation fields D, H in the usual 
Maxwell electrodynamics. The constitutive equations have the form 

antisymmetric in Zi/ 2 as well as qi 9253 (the indices of k have been lowered with the gen- 
eralized Minkowski metric tensor; we shall treat other tensors in the same way in the 
following). It is useful at this point to distinguish between the space-time elements 
and the elements f^. Let us assume that the tensor introduced in Eq. (24) does not mix 
these fields (for n M5 , if £ t i5 qiq2 q 3 has qiq2Q3 = Oifi'j, then the components of k that enter 
are of the form k a p 1 = Sap-y^sf^ 5 only; similarly, for n^, only the components £^ va f3b 
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can occur, and /c a/35 connects only to the components f\ a of the field tensor). As we have 
pointed out above, the vector is physically distinguished from the antisymmetric tensor 
f^v in the inhomogensous field equations, since the source terms break the higher symme- 
try of the homogeneous field equations. The assumption that the constitutive equations 
do not couple these components results in a simpler system to analyze, although (as for 
Hall type effects in the non-relativistic theory) it is conceivable that the more general case 
could occur. 

We introduce the new set of fields: 

= \e^ vXa h^ (25) 

so that 

n \a = 2S Xaa ' 35 b a p. (26) 

On the zero mode, the fields correspond to the dual Maxwell fields; in this theory they 
play a role analogous to the B fields in the Maxwell theory. In a similar way, the / M 5 fields 
are analogous to E. The part of the tensor Ei x i 2 qi q2 q3 connecting the fx5 fields is discussed 
below. 

Working with these fields enables us to construct the equations in a form which, as 
we shall show, generalizes the Maxwell theory to a form where the invariant time r plays 
the role of t and spacetime plays the role of space. This analogy helps to interpret the 
physics and it distinguishes between the familiar physical quantites and the new fields 
/ M 5. Substituting these fields in (21) and (22), we find 

d 5 b^ + d a e^ x f X5 = 0, (27) 



analogous to (Id), and 

d a V" a = 0, (28) 
analogous to (lc). For the spacetime excitation fields we define 



h - I > 



and we get from (23), for Z 2 = /z, 



d 5 n^ - l - e ^d a h Xv = (29) 

analogous to (16), where we have used 

e af3vS s^ = -2(5Z5^-5^). (30) 
To complete the set of equations, we note that for li = 5, we get from (23) 

dan^ = -f, (31) 
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analogous to (la). 

To obtain a mass-energy conservation law for the fields, we multiply (29) by and 
(27) by h^, and then combine them, obtaining 

[U 5 d 5 n^ + l -h^d T b» v ] + l -e^ Xv d a {f^h Xv ) = -fj^. (32) 

Assuming the dielectric tensor reduced into the ju5 and [iv subspaces is symmetric (the 
relations of n^ 5 to / M 5 and to go by the contraction £e; the exclusive property of 
indices of e then imply simple conditions on £ for the symmetry of these forms) we can 
write (32) as 

\ds [f^ 5 + \h^] + \e^dM^hx,) = -J M ^5- (33) 
From the Stueckelberg Hamiltonian 17 

K = 7^(P M - e a M )(p^ - e a M ) - e a 5 , (34) 

where eo differs from the usual electric charge by a dimensional constant, one can derive 
the relativistic Lorentz force 13 ' 17 ' 21 ' 24 

M^ = e r^ + e r 5 , (35) 

from which it follows that 

Mi„i" = M-^-(i^) = xj^ (36) 
ar 

Since, in the Stueckelberg theory, the Hamilton equations imply that 

and hence the r derivative in the central equality of (36) corresponds to a change in the 
mass-squared {j>^ — eoa M )(p M — eoa M ) of the particle. It therefore follows that j^f^s, is the 
rate of change of the mass of the particle. We therefore identify s a = \s a ^ Xu f ^h\ u as the 
analogue of the Maxwell Poynting vector. This Poynting 4-vector is the mass radiation 
of the field. We see, furthermore, that \ [f^n^ 5 + \h^ u b^ ul \ is the scalar mass density of 
the field (its four integral is the dynamical generator of evolution of the non-interacting 
field 17 ' 19 ). 

We now introduce the eikonal approximation, i.e., set 

fhl 2 ( x , T ) = fhl 2 ( x ) e *P iK ( T ~ *<») 

for large k. In the absence of sources the 5D-Maxwell equations (27), (28), (29), (31) take 
the form (for large k) 

^-£^V/A5 = 0, (37) 
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VJT" = 0, (38) 
M 5 + }. e ^ pahxv = o, (39) 

Pan 5a = 0, (40) 



" ' 2 



where p a = d a ~^> . 

We now relate the direction of p^ to the polarization of the fields. We write the "cross 
product" of n and b (analogous to the cross product of D and B in Maxwell's theory): 



e^xn vh b aX = 2n» 5 f„ 5 p^ (41) 



or 

Vil = 2n^f a5 S ^ XnU5baX ' 

where we have used (37) and (40). It is clear that since p M and the Poynting four- vector 
are cross products of tensors which are not necessarily aligned in the same four-directions, 
they are in general not parallel to each other (in space-time) due to the anisotropy of the 
medium, i.e., the wave normal and radiation flow directions are not, in general, the same. 

The relations (37) — (40), along with the constitutive relations relating n Ml/ , f a \, and 
n^ 5 , /o- 5, provide relations analogous to (6) characterizing the possible field strengths of 
the eikonal approximation in terms of properties of the medium. We shall not treat these 
relations here, but discuss the mass-radiation flows, along the rays, on spacetime geodesies 
in the interesting special case where W v = b^ u , which is analogous to the case of mate- 
rials with \x — 1 in Maxwell's electromagnetism. This case is interesting since, although 
the space is empty in the usual sense (in analogy to E = D , B = H) , the dielectric ef- 
fect involving the / M 5 components can drive the radiation on curved trajectories, i.e., the 
corresponding spacetime can have a non-trivial metric structure. 

We multiply (37) by \s al3 » v pp. We then use (39) and (30) to obtain 

nj - p a pPf/35 + P^faS = (42) 

Defining the reduced dielectric tensor S a ^ as the part of the general dielectric tensor 
which connects only the a5 components of the fields, i.e., 

n a5 = £Jfps, (43) 

the condition (40) then implies that £q//35 cannot be in the direction of p a (unless it is 
lightlike). We obtain from Eq. (42) 

(S a 13 - PaP^J + p a p p ) fps = 0, (44) 

where we have chosen the negative sign for the signature of the fifth index, n a ^ = — n a 5 
[with this choice the flat space limit, for which £ a @ — 5 a " , Eq. (44), with (40), admits only 
spacelike p a ; for positive signature of the fifth index, in this limit, p a would be timelike]. 
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Eq. (41) has a solution only if the determinant of the coefficients vanishes (a similar 
calculation in which the field strengths f^ v enter in place of results in the same condition 
on these coefficients, as it must). It is somewhat simpler to work with the eigenvalue 
equation (44) . Assuming as before that this dielectric tensor is symmetric, we can work in 
a Lorentz frame in which it is diagonal. In this frame we have (for the transformed fields) 

Note that in the isotropic case for which all of the E a are equal, one obtains ppf 135 = 0, 
and the metric becomes conformal, i.e., one obtains the condition 

= -i. 

where rf" v is the flat space Minkowski metric (—1, 1, 1, 1). 

Multiplying the equation (45) on both sides by p a , and summing over a, one obtains 
the condition (p 2 = p^p^), 

= K=y^ + 7 A^ + 7 A^-y^ + l. (46) 

Ei -p 2 E 2 - p 2 E 3 -p 2 E - p 2 

This condition determines, in this case, the Fresnel surface of the wave fronts. 
It then follows that 

dK 2p» ^8K_ 

dp„ ~ E^-p 2 + P dp 2 ' { ' 

Calculating the scalar product of (45) and (47) one then obtains 
dK _ 

(f) 2 Po 2 OK ' ( 48 ) 



-2 



J \ = ^ 3 (^-P 2 ) 2 (E°-p 2 ) 2 dp 21 



Multiplying the expression (37) for b^ft^) by (47), the contribution of the second term 
of (47) vanishes since the Levi-Civita tensor is antisymmetric; the first term, according to 
(45), is proportional to / a5 , and vanishes for the same reason. It therefore follows that 

OK 

t—W = 0. (49) 
dpv v ' 

Since the scalar product of J^- with both W and f^ 5 is zero, it is proportional to their 
"cross product" i.e., it is parallel to the Poynting vector. To make the proof explicit, it is 
convenient to define = J^-, Hi = —h j, F t = f i5 , and D l = e^ k hjk (the space index 

may be raised or lowered without changing sign in our Minkowski metric). In this case, 
the conditions V^h^ = and V^f^ become 

V°H - V x D = 

-F°/ 05 + V-f = (50) 
V H = 0, 
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where we have used boldface to represent the space components of the vector. In these 
terms, the Poynting vector is given by 



S° = D • f 
S = / 05 D + f x H 



... TT (51) 



Taking the cross product of f with the first of (50), one obtains 

V°(f x H) = V(f • D) - D(f • V). 

For V° 7^ 0, one may substitute this into the second of (51). The / 05 D term, with the 
help of the second of (50), cancels, and we are left with 

S = — V. 

v° 



v 

case V° = 0, the second and third of (50) imply that 



It then follows that = t^V^, i.e., is proportional to the Poynting vector. For the 



V f = V H = 0, 



i.e., if V 7^ (the case = is exceptional in the eikonal approximation), it must be 
proportional to f x H. From the first of (50), we see that V x D = 0, and if D ^ 0, it must 
be proportional to V. The space part of S^, from the second of (51) is then proportional to 
V. Under these conditions, the time part of S^ 1 vanishes, and therefore we again obtain the 
result that is proportional to 5 M . If D = 0, then S° = and, since V is proportional 
to f x H, it again follows that is proportional to S^. 

From this point, one may follow the same procedure used in the case of Maxwell's 
electromagnetism (Eqs.(14) to (18)) to obtain the Hamiltonian flow corresponding to the 
admissible modes. The Lagrangian associated with the Hamiltonian (46) satisfies a mini- 
mal principle, from which it follows that the Hamiltonian flow is geodesic on this manifold. 
Replacing f 135 in (44) by (^" 1 ) /3 7 n 75 , one obtains 

{5% - M% M; yy> 75 = 0, (52) 

where 

M a 7 „ " = 8 a 7 - ^ "(S-Y 7 - (53) 

The condition (40) implies that the solutions can lie only in the hyperplane orthogonal to 
p 17 . The projection of the matrix M^p^p u (in the indices a, 7) into the three dimensional 
subspace orthogonal to p a is symmetric, and can therefore be diagonalized by an orthog- 
onal (or pseudo-orthogonal) transformation in three dimensions. In fact, the Gauss law 
and a gauge condition restrict the polarization degrees of freedom to three 19 (the eikonal 
approximation is far from the zero mode, which corresponds to the Maxwell limit, for 
which only two polarizations survive) and hence one finds three geodesies. 
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For p a timelike, one can choose a Lorentz frame in which the eigenvalue condition 
(52) has the form 

{5%-M- (p'°) 2 }^ 5 = 0, (54) 

and for p a spacelike, 

{5% - M^ 3 (p" 3 ) 2 }n"^ = 0. (55) 

In each of these cases, the matrix can be diagonalized under the little group acting in the 
space orthogonal to p a (leaving it invariant). For the lightlike case, up to a rotation, p a has 
the form (p, 0, 0, p). The remaining matrix may then be diagonalized under SO(2) rotations, 
to obtain just two geodesies, corresponding to the polarization states of a massless Maxwell- 
like theory. This special limiting case will be investigated in detail elsewhere. 

With the same procedure as applied to the Maxwell case treated above, with H 
replaced by K, and the space indices replaced by spacetime indices, one finds the symplectic 
structure of the flow of matter in space time. 

It has been shown by Kline and Kay 2 , as discussed above, that for the three dimen- 
sional Maxwell case, the Hamilton equations resulting from the eikonal coincide with the 
geodesic flow generated by the resulting metric (recall that the direction of momentum 
associated with both eigenstates is the same); a similar proof can be applied to the 4D 
cases we have studied here. 

4. Eikonal of the Relativistic Stueckelberg-Schrodinger Equation 

It is interesting to apply the eikonal method to a relativistic Schrodinger equation in 
a medium which is not isotropic, for example, in a crystal with shear forces 25 , with locally 
varying band structure (as in a crystal under nonuniform stress, or near the boundaries 
or impurities). In this case, the rays are directly associated with the (probability) flow of 
particles. The eikonal eigenvalue condition is one dimensional in this case, since the field 
is scalar. For an analog of this structure (corresponding, for example, to a distribution 
of events in a potential periodic in both space and time) in four dimensions described 
by a relativistically covariant equation of Stueckelberg-Schrodinger type, the metric one 
obtains is a spacetime metric, and the geodesic flow is that of the quantum probability 
for the spacetime events (matter) described by the Stueckelberg wave function. We show 
elsewhere 26 , moreover, that such an equation may be derived from a relativistic generaliza- 
tion of a stochastic procedure analogous to that of Nelson 27 , with non-trivial correlations 
between the stochastic motions in different directions of spacetime. This simplest analog 
of the nonrelativistic problem is given by 

i^ T [x) = d lx E lJtV d v ^ T {x), (56) 

where corresponds to the effect of the medium, and is assumed to be symmetric. The 
Schrodinger current is then 

3r(x) v = -itySS^d^r ~ VV<V<^VV*). (57) 

In the eikonal approximation, for which the frequency associated with r (essentially the 
total mass of the particle 15 ) is large, one obtains the condition 

K = £^p u -1 = 0, (58) 
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analogous to the Fresnel surface condition (46) for the optical case. It is clear that dK/dp^ 
is in the direction of j£. This implies that K is the operator of evolution for the dynamical 
flow of particles, corresponding to the rays. It follows from the Hamilton equations that 
the flow is geodesic, where £^ v is the metric. 

5. Summary and Discussion 

We have reviewed the eikonal treatment of Maxwell waves in an anistropic medium in a 
form which exhibits the structure of a canonical mechanics on a three dimensional pseudo- 
Riemannian manifold, where there are two metrics, one for each admissible polarization 
mode. The canonical momenta for each of the corresponding geodesies have the same 
direction, but differ in magnitude. 

We then studied the eikonal approximation for the 5D wave equation for the gener- 
alized Maxwell fields implied by the requirements of gauge invariance of the Stueckelberg 
manifestly covariant quantum theory. One finds in this case the structure of a canonical 
mechanics on a four dimensional pseudo-Riemannian manifold, where there are in general 
three metrics for each of the timelike or spacelike possibilities for the canonical momenta, 
and two for the limiting lightlike case. No motion of the medium is required. This canon- 
ical mechanics has the form of Stueckelberg's covariant mechanics, where the form of the 
Hamiltonian describes the Fresnel surfaces, and the evolution generated by the Hamilton 
equations carries mass flow along the geodesies. 

We finally studied briefly the eikonal approximation to an non-isotropic Stueckelberg- 
Schrodinger equation written in a form analogous to the non-relativistic theory in a crystal 
as it appears in a Brillouin zone, where the local energy surface carries curvature (an exam- 
ple of a model providing an analogy to a crystal in spacetime is that of an electromagnetic 
field in a resonant cavity), or, as derived from a relativistic Nelson type stochastic pro- 
cedure with non-trivial correlations 26 . The eikonal approximation corresponds to (total 
system) mass large compared to typical Hamilton-Jacobi momenta. In this case the Hamil- 
tonian flow is in the direction of the eikonal approximation to the quantum mechanical 
current, so that the analog of the Fresnel surface condition is the generator of evolution 
for the flow of particles, corresponding to the rays. The Hamiltonian flow is geodesic with 
respect to the pseudo-Riemannian metric identified with the effective mass matrix. This 
construction provides a model for an underlying quantum mechanical structure for classical 
dynamical motion along geodesies in a pseudo-Riemannian manifold, an analog to general 
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